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ABSTRACT 


An iodine-operated 200-W Hall thruster plume has been simulated using a hybrid-PIC 
model to predict the spacecraft surface-plume interaction for spacecraft integration purposes. For 
validation of the model, the plasma potential, electron temperature, ion current flux, andion number 
density of xenon propellant were compared with available measurement data at the nominal 
operating condition. To simulate iodine plasma, various collision cross sections were found and 
used in the model. While time-varying atomic iodine species (i.e., |, I*, ?*) information is provided 
by HPHall simulation at the discharge channel exit, the molecular iodine species (i.e., lz, l2*) are 
introduced as Maxwellian particles at the channel exit. Simulation results show that xenon and 
iodine plasma plumes appear to be very similar under the assumptions of the model. Assuming a 
sticking coefficient of unity, iodine deposition rate is estimated. 


INTRODUCTION 


Hall-effect thrusters (HET) are relatively low-cost electric propulsion devices that can be 
utilized for a wide range of ambitious long-duration NASA missions. Xenon (Xe) is an attractive 
propellant because it does not chemically react with spacecraft surfaces, and has larger mass and 
ionization cross section compared to other inert gases, yielding higher thrust for a given input power 
(T/P). However, Xe is relatively expensive and must be stored inlarge high-pressure tanks [1]. Due 
to these limitations, a variety of alternatives have been proposed, and iodine has been identified 
as an attractive alternative propellant for HETs, especially for small satellites (<3800kg) and 
exploration class missions [2]. lodine is almost as heaw as, is easier to ionize than, and has about 
three times the storage density of Xe [3]. The feasibility of a Hall thruster operation using iodine 
propellant has been proven through experimental testing of Busek’s 200-W Hall thruster (BHT- 
200), which yielded stability and performance comparable to that observed with Xe, higher thrust 
to power ratio at higher power operating condition, and lower beam divergence at certain operating 
conditions [8, 4]. 


Despite its advantages, using iodine as a propellant has one critical risk, which is 
associated with its high reactivity. Reactive atomic iodine (I) may potentially cause contamination 
or corrosion of spacecraft surfaces and solar arrays, which is a concern for spacecraft system 
integration. Since the majority of Hall thruster experiments that utilize iodine have focused on the 
thruster and subsystem performance, there is adearth of detailed knowledge of physical processes 
occurring in the plume of the thruster and thus their effect on spacecraft surfaces. Moreover, due 
to the contribution from molecular iodine (lz), many different types of reactions may occur in some 
regions—including dissociative ionization, electron attachment, and inelastic energy exchange— 
that are very different from the nature of Hall thrusters operated on Xe. Although there have been 
extensive experimental efforts in using iodine as a Hall thruster propellant, results from modeling 
work associated with the iodine thruster plume and its effect on spacecraft surfaces have not been 
published to date. Modeling the iodine plasma plume andits interaction with the spacecraft surfaces 
is necessary in order to further understand its use as apropellant and mitigate any potential harmtul 
surface-plume interactions. 


This study is concerned with numerical simulations of the iodine plasma plume generated 
by BHT-200 Hall thruster (shown in the left of Figure 1), which would be the first iodine-operated 
Hall thruster that will be flown in space in the lodine Satellite ((SAT). The basic configuration of 
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ISAT, including the BHT-200 thruster and other components of the propulsion system are shown 
in the right of Figure 1, and more details about the iSAT mission canbe found in Ref. [1]. Since the 
iodine plasma properties of the BHT-200 thruster plume have not been measured, the simulation 
is performed first using xenon as a propellant for a validation of the model. To simulate the iodine 
plasma plume, various collision cross sections were found and used in the hybrid PIC model, 
Draco. In a Hall thruster plume, iodine exists as both atomic (i.e., |, I*, I?*) and molecular (i.e., Iz, 
l2*) species. Due to these mixtures, complex physical processes may occur, and thus modeling an 
iodine plume can be challenging and requires some assumptions, which are discussed in the next 
section. Finally, the iodine particle flux on the spacecraft surface is found to estimate its deposition 
rate. 


Figure 1. Left: Busek’s 200 W Hall thruster (BHT-200). Right: Basic configuration of iSAT including the 
BHT-200 Hall thruster in the front and the solar arrays in the back of the satellite. 


MODEL DESCRIPTIONS 


The model employed in this study, Draco [5], is a software framework developed under 
COLISEUM framework by Air Force Research Laboratory (AFRL). Draco is a 3-D parallel hybrid 
code that was developed to study electric thruster plumes and their interaction with spacecraft 
components. In Draco, the surface is defined using a triangular mesh, while a Cartesian volume 
mesh is used for the Draco’s computation domain. For an illustration of the surface and volume 
meshes used in this simulation, Figure 2 shows 1) the surface mesh of the iSAT spacecraft and 
the BHT-200 Hall thruster, created using Cubit, and 2) the rectilinear volume mesh generated by 
Volcar—Draco’s built-in mesh generator—with a contour plot of the number density of Xe*, on the 
left and right sides, respectively. 
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Figure 2. The geometry and surface mesh of the iSAT spacecraft and BHT-200 thruster, created using 
Cubit (left), and the cross-section view of the volume mesh with ion number density contour (right). 


PLASMA DYNAMICS 


In the current model, the transport of heaw species are simulated using the particle-in-cell 
method (PIC) [6] combined with either direct simulation Monte Carlo (DSMC) [7] or Monte Carlo 
Collision (MCC). The motions of charged particles in plasmas are governed by the Lorentz force F: 


F=q(E+vxXB) (1) 


where E is electric field, B is magneticfield, and v is particle velocity. In Hall thrusters, the magnetic 
field is not strong enough to affect the motion of ions. 


In the PIC method, the charge density at each grid point is calculated from the distribution 
of particles in each cell. From the charge density, quasi-neutrality (n, ~ n,;) is assumed, and the 
plasma potential is obtained by solving the Boltzmann relation that can be derived from electron 
momentum conservation equation assuming thermal equilibrium, constant electron temperature, 
and no magnetic fields: 


KgTer Ne 
es aera on (2) 
where ¢ is electrostatic potential, T, is electron temperature, n, is electron number density, and 
subscript r represents a reference value. Instead of assuming a constant temperature in the entire 
domain, a polytropic temperature model can be used to provide more accurate representation of 
the electron temperature field, which results in the following modified Boltzmann relation: 
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where y is specific heat ratio. The electric potential is used then to calculate the electrostatic field 
and force on the particle from E = —V®@. Then, the particle positions and velocities are updated by 
integrating Eq. (1) at each time step. 


COLLISION DYNAMICS 


In a Hall thruster plume, the two most important collision mechanisms are the momentum- 
exchange (MEX) and the charge-exchange (CEX) collisions between the neutral and ion species. 
The current study employs the MCC method [6] to simulate both MEX and CEX collisions between 
neutral-ion pairs, and the variable hard-sphere (VHS) model [7] to simulate the MEX collision 
between neutral-neutral pairs. In the MCC method, a source particle collides with a target “cloud,” 
instead of an individual target particle. A collision probability is given by Birdsall [6] as: 


P =1-—exp(—n, cgAt) (4) 


where n,, is the number density of target particles, o is the collision cross section, and g is the 
relative velocity. If this probability is larger than a random number, a collision occurs and the type 
of collision is determined. 


If aCEX collision occurs between an ion and a neutral cloud, the velocities of the collided 


ion are replaced with random values obtained by sampling the Maxwellian distribution at the neutral 
temperature. The Maxwellian distribution is given as: 


fu = pe (Ye es 5) an (5) 


where k, is the Boltzmann constant, T is temperature, m is particle mass, R; is arandom number 
between 0 and 1, and M is the number of R;. 


To account for the strong forward-scattering tendency of fast ion and slow neutral 
collisions, scattering angles are calculated using adifferential cross section. A fit to an experimental 
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CEX differential cross section measurement for xenon in a Hall thruster [11] is found and 
implemented by Scharfe [8]. 


CROSS SECTION MODELS 


As shown in Eq. (4), a collision cross section is necessary to calculate a collision 
probability. For neutral-neutral collisions, the VHS cross section [7] is employed: 


Pio... 3 
Oyns(Xe,Xe) = gaa m (6) 


20 
where g is the relative velocity and w = 0.12 is related to the viscosity temperature exponent for 
xenon. For modeling CEX collisions, Hall thruster modeling community often uses the following 
semi-empirical form of cross sections: 


Ocex = A— B log(E) (7) 


where A and B are fitting parameters to experimental data and E is the collision energy in eV in the 
laboratory frame. For xenon collisions, the CEX collision cross sections between a Xe neutral and 
single- and double-charged Xe ions were measured by Miller ef. a/. [9], and the fitting coefficients 
to Miller's experimental data using Eq. (7) are summarized in Table 1. 


Table 1. Coefficients in Eq. (7) for xenon and iodine CEX collision cross sections. 


A B 
Xe-Xe* 87.3 13.6 
Xe-Xe*t 45.7 8.9 


Modeling an iodine plasma plume can be more complicated than xenon because there are 
more reaction and collision types to account for due to the presence of molecular iodine species, 
i.e., 1, and IZ. The likelihood of various collisions occurring depends on their collision cross 
sections. The present work does not attempt to solve all rate equations to determine which 
collisions are important to consider. Instead, just as in a xenon-operated Hall thruster plume, it is 
reasonable to assume that an efficient long-range CEX collision is important also in an iodine- 
operated Hall thruster plume, because the CEXcollision produces slow ions and fast neutrals that 
affect the thrust performance and lifetime. Since iodine exists as in both atomic and molecular 
states in a Hall thruster plume, we need to consider the CEXcollisions between all of the following 
cases: 1) atomic neutral and atomic ion (I-*), 2) molecular neutral and atomic ion (/,-/*), and 3) 
molecular neutral and molecular ion (/,-/7). Other species combination could be considered for 
more accurate collision modeling, but only these three types of collisions are considered in the 
current study. Recently, the CEX collision cross sections for /,-I1* and I,-17 were measured by 
Hause et. al. [10]. The fitting coefficients to Eq. (7) are found for I,-1+, as A= 66 Aand B = 4.7 A, 
as summarized in Table 1. Since a fit for I,-17 data to Eq. (7) was not determined in Ref. [10], we 
determined an analytical fit to this CEX cross section data using the following form: 


Ocex(*,1,) = c,log3(E) +c, log?(E) + cz log '(E) + cy (8) 
where c,,Cz,C3, and c, are fitting coefficients defined as in Table 2 below. 
Table 2. Values of the fitting coefficients in Eq. (8). 
C4 C2 C3 C4 
-0.47 3.5 -9.0 82.0 


Since a measurement of the cross section of J-1* is not available, a semi-empirical formula 


that relates the cross section with its ionization potential, ¢,, and impact velocity, v, by Sakabe and 
Izawa [11] is used to calculate the unknown cross section: 


-1.5 
o(v) = [A-B log,)(v)] (=) (9) 
Io 

where A = 1.81 x 10°“, B = 2.12 x 10-*°, and e,, = 13.6 eV is the ionization potential of hydrogen. 
For iodine, the first ionization potential of atomic iodine was used (i.e., €, = 10.5 eV). In order to 
verify this cross section model, the cross section of xenon is calculated using Eq. (9) and is 
compared with Miller’s cross section for xenon as shown in Figure 3. Sakabe’s formula predicts 
Miller’s cross section for xenon very well. Thus, the cross section for J- /* is found using Sakabe's 
formula and employedin current study. AllCEXcross sections thatare modeled in Draco are shown 
in Figure 4. The CEX cross section for J- 1* is lower than Xe- Xe* for increasing impact energy, so 
fewer collisions are expected for /- 1* than Xe- Xe*. On the other hand, /,- 1/7 is higher than Xe- 
Xe* for increasing impact energy. 
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Figure 3. Charge-exchange collision cross section as a function of impact energy. Sakabe’s formula 
predicts Miller’s semi-empirical fit very well. 
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Figure 4. Charge-exchange collision cross sections for all species modeled. 


KEY ASSUMPTIONS USED INMODELING IODINE PLASMA PLUME 


Performing a plasma plume simulation using Draco requires neutral and ion particle 
information at the thruster exit as an input condition. The effect of input source types in Draco have 
been studied before [12, 13]. Among other source options, utilizing HPHall Spawn to provide the 
particle information to the plume was proven to be the most effective option. Thus, we 
simultaneously run HPHall [14] to model the Hall thruster discharge plasma and to provide an 
accurate time-dependent particle information to Draco at each time step. Figure 5 shows the 
simulation domain of HPHall, with a contour plot of an electron number density. HPHall sends the 
particle information to Draco along a user-defined lambda line (not shown). 
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Figure 5. Contour plot of electron number density from HPHall simulation, ran simultaneously to 
provide detailed input to Draco. 


Due to the mixture of atomic and molecular iodine species, a number of reactions involving 
lz may occur in the thruster and the plume, including dissociation, ionization, electron attachment, 
inelastic energy exchanges (e.g., rotational and vibrational), and recombination. Especially, since 
lz has alow dissociation energy threshold (1.54 eV) and alarge ionization cross section (peak cross 
section = 12.3 A’), a dissociative ionization process is likely to occur in the ionization region of the 
discharge channel. This process should be modeled in HPHall for a physically accurate model of 
the iodine plume. Since HPHall currently models only the neutral atom, singly - and doubly-charged 
ions, implementing additional reactions in HP Hall involves modifying significant portion of HP Hall. 
More importantly, a detailed reaction modeling is not the focus for this work. Rather, the goal is to 
provide a first-order approximation of the iodine particle flux on spacecraft surfaces using the 
numerical tools available to us at this stage. 


For the iodine plume simulation, only J, /*, and 12+ species are simulated in HPHall, and 
I, and I} are injected as separate inflow species at the thruster exit plane assuming Maxwellian 
velocity distributions, resulting in all five iodine species (/, /*, I?*, I,, and 1+) handled by Draco in 
the plume. The mass flow rate of Iz* has been estimated based on measured mole fraction data 
shown in Table 3, and the mass flow rate of lz is assumed to be 10% of the total neutral flow rate, 
which is calculated from the mass utilization efficiency of 85.3% [3]. 


All simulations are performed at the nominal operating condition of the BHT-200 thruster, 
as shown in Table 4, with the background neutral number density equivalent to the vacuum 
chamber pressure of 5x10 Torr. The MEX and CEX collisions are ignored in HPHall simulations 
but included in the plume simulation using Draco. The mass and ionization cross section for each 
species was appropriately incorporated in HPHall. 


Table 3. Xenon and iodine ion species distributions in the BHT-200 Hall thruster plume [3]. 


Species Mole Fraction, Xe Mole Fraction, | 
It 0.029 

Xe* I* 0.975 0.953 

Xe?* 12+ 0.021 0.015 

Xe** 13+ 0.004 0.003 


Table 4. Summary of nominal operating conditions for the BHT-200 Hall thruster. 


Xenon [12] lodine [3] 


Discharge voltage (V) 250 250 
Discharge current (A) 0.75 0.74 

Anode mass flow rate (mg/s) 0.84 0.82 
Cathode mass flow rate (mg/s) 0.098 0.096 
Mass (propellant) utilization efficiency 0.981 0.853 


RESULTS AND DISCUSSION 


As a model validation, the first comparison is made between the simulation results and the 
laboratory measurements of xenon plume plasma parameters of the BHT-200 thruster from Ref. 
[12]. The left side of Figure 6 shows angular profiles of the plasma potential at a radius of 10 cm 
away from the thruster exit plane, compared with experimental data. The simulated plasma 
potential value shows an excellent agreement with the experimental data, where both results 
monotonically decrease from 19 V to 2.0 V, from 0 degree to 90 degrees. The right side of Figure 
6 compares the angular profiles of the simulated and measured electron temperature at 10 cm 
away from the thruster exit plane. The electron temperature predicted by the model decreases from 
3.8 eV to 0.94 eV from 0 to 90 degrees, while the measured data decreases from 2.6 eV to 0.36 
eV. Despite the over-estimation of Te by 1 eV, the overall trend agrees well with the measured 
data. The small discrepancy in the magnitude of the electron temperature is likely due to the 
limitation of the polytropic temperature law and/or the choice of the reference values for the 
Boltzmann relation. 
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Figure 6. Angular profile of the plasma potential (left) and electron temperature (right) at a radius of 
10 cm away from the centerline of the thruster exit plane. 


Draco 
Oo Measurement 


Draco 
oO Measurement 


£ Pn 
< “10° 
= : E 
= eo5 r= 
5 oO 77 
a1 ale S 
o 15 
= oO 10 
o [Oo 3 
+ — 3 
e On o E 
= 
8 Do 
10" 
102 gs is 
20 40 60 80 10 
Angles (deg) 0 20 40 60 80 
Angles (deg) 


Figure 7. Angular profile of the ion current density from an RPA probe (left) and the ion number 
density (right) ata radius of 10 cm away from the centerline of the thruster exit plane. 


The left side of Figure 7 shows the numerical angular RPA trace, compared with the 
collected ion current density by an RPA probe. The ion current density is over-estimated near the 
thruster centerline (<5 deg.) and is under-estimated at larger angles (> 20 deg.), especially around 
60-70 degrees where the charge-exchange wing is observed from the measurement. The 
comparison shows that the model predicts a more collimated plume than the measurement at this 
operating condition. The right side of Figure 7 compares the ion number density predicted by the 
model with experimental data. While the model provides reasonably good magnitudes in general, 
the overall shape is not well captured. Although the discrepancies in ion properties could be 
minimized by using better choices of the reference values for the Boltzmann relation in Draco and 
anomalous mobility parameters in HPHall, it was decided that the model provides a reasonable 
agreement with experimental data overall for the purpose of this work, which is to estimate the flux 
of iodine on the spacecraft surfaces. 


Next, we compare the xenon plume simulation results with the iodine. In Figure 8, number 
densities of the total neutral and ion species are compared on the left and the right, respectively. 
From the qualitative comparison, both Xe and! showsimilar contours ofthe neutral and ion number 
densities overall. Xe and | have similar mass and ionization potential, and there are very few 
molecular iodine species in the plume. The sum of the molecular species (not shown) have at least 
an order of magnitude smaller current and number densities than the singly-charged atomic iodine 
ion (I*). Similar observation is made in the plasma potential and ion current density for Xe and | in 
Figure 9. These results are not surprising because similar performance and plume efficiencies 
between Xe and! were measured at this operating condition before [4]. Moreover, the mass flow 
rates of the molecular species were based on experimental data. 


The contour plots of the total neutral and ion particle fluxes on the surface of the iSAT 
satellite are shown in Figure 10. The fluxes of both neutrals (I and |2) and the ions (I*, I?* and le*) 
on the surface of the iSAT decrease away from the thruster exit channel because the neutral 
number density in the plume decreases as 1/r?, where r is the distance from the thruster exit. 
However, the fluxes appear to be slightly higher towards the outer edges of the front surface of the 
iSAT body than the inner region. Moreover, the front of the solar array shows the most amount of 
particle fluxes. The highest neutral andion fluxes on the entire spacecraft are observed on the front 
pole pieces of the thruster. Due to the low density of neutral particles, the surface flux contour 
appears to be coarse, which would improve with increased run-time. 
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Figure 8. Comparison of total neutral (left) and ion (right) number densities of xenon (Xe) and iodine 
(l) species. 
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Figure 9. Comparison of the plasma potential (left) and total ion current density (right) of xenon (Xe) 
and iodine (I) species. 


Based on these results, the highest neutral and ion particle fluxes on the solar array are 
calculated as approximately 3.5x10' m’s"' and 1.0x10'® m*s"', respectively. Assuming 100% of 
the total incident particles sticks to the surface of the solar array and a constant rate of deposition, 
whichis an extremely high conservative estimate, the deposition rate ofiodine particles onthe solar 
array would be approximately 0.007 um/hr, which translates to the deposition per unit of 
approximately 0.34 mg/cm? over the entire thruster operation duration (~10 hrs). 


However, not all iodine particles hitting the solar panel surface will deposit or chemically 
react to the surface, given high vapor pressure of iodine at the temperature of the solar panel. The 
portion of particles that actually react to the surface will depend on the surface properties of the 
solar panel, which is beyond the scope of this work. Moreover, the incident iodine flux on the 
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thruster surface is almost an order of magnitude higher than the flux on the spacecraft body or solar 
array, which could be a concern in the future for much longer thruster operation duration than the 
iSAT mission. 
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Figure 10. Contours of the total neutral (left) and ion (right) particle fluxes of iodine on the iSAT 
surface in particles/(m?s). 


SUMMARY AND CONCLUSIONS 


In this study, both xenon- and iodine-operated BHT-200 Hall thruster plumes were 
simulated using Draco to estimate the incident particle flux on the spacecraft surface of the iSAT 
satellite. The model was first validated using traditional Xe propellant by comparing with available 
data. Then, the iodine propellant—including /, /*, 1?*, I,, and I+ species—was simulated using the 
collision cross sections found from literature. The atomic iodine species (I, 1*, and /?*) were 
simulated using the HPHall to pass the particle information from the discharge channel to the 
plume, while the molecular species (/,, and I7) were introduced at the discharge channel exit 
assuming Maxwellian velocity distribution functions. Under these assumptions, the iodine 
simulation results looked very similar to xenon, likely because of the similar mass and ionization 
cross section, and low molecular species density in the plume. 


Not all iodine particles hitting the solar panel surface will deposit or chemically react to the 
surface, given high vapor pressure of iodine at the temperature of the solar panel. However, the 
sticking coefficient was assumed to be unity as a highly conservative estimate. Based on this 
assumption, the highest total iodine flux on the solar array was calculated as approximately 
4.5x10'® m*s', which means the deposition per unit area would be approximately 0.34 mg/cm* 
over the entire thruster operation duration. In reality, only some portion of iodine colliding with the 
surface may chemically react with the surface, which could reduce the efficiency and/or 
performance of the solar panel. How many particles actually react to or reflect off the surface will 
depend onthe surface properties ofthe solar panel, whichis beyond the scope ofthis work. Finally, 
for more physically accurate simulation of iodine plasma plume, one needs to model the detailed 
oe discussed in this paper, especially the dissociative ionization, inside the discharge 
channel. 
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